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Anomalies of the infrared-active phonons in underdoped YBCO 
as an evidence for the intra-bilayer Josephson effect 
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The spectra of the far-infrared c-axis conductivity of underdoped YBCO crystals exhibit dramatic 
changes of some of the phonon peaks when going from the normal to the superconducting state. We 
show that the most striking of these anomalies can be naturally explained by changes of the local 
fields acting on the ions arising from the onset of inter- and intra-bilayer Josephson effects. 

PACS Numbers: 74.25.Gz, 74.72.Bk, 74.25.Kc, 74.50. +r 
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The essential structural elements of the high-Tc su- 
perconductors are the copper-oxygen planes which host 
the superconducting condensate. Many experiments, and 
also some theoretical considerations, suggest that these 
planes are only weakly (Josephson) coupled along the c- 
direction. Studies of the c-axis transport |jl|,H and those 
of the microwave absorption and the far- infrared 

c-axis conductivity revealing Josephson plasma res- 
onances, have established that Josephson coupling indeed 
takes place for planes (or pairs of planes) separated by in- 
sulating layers wider than the in-plane lattice constant. 
It is not fully understood why the coupling is so weak 
and it is debated whether this is related to the uncon- 
ventional ground state of the electronic system of the 
planes causing a charge confinement and/or to the 
properties of the insulating layers. In this context, it is of 
interest to ascertain whether the closely-spaced copper- 
oxygen planes of the so-called bilayer compounds, like 
YBa2Cu30j,, are also weakly (Josephson) coupled. 

In this paper we show that the far-infrared spectra of 
the c-axis conductivity of underdoped YBa2Cu30j, with 
6.4 < y < 6.8 may provide a key for resolving this inter- 
esting issue. The spectra exhibit, beside a spectral gap 
that shows up already at temperatures much higher than 
Tc Ipj-ll^, two pronounced anomalous features 
Firstly, at low temperatures a new broad absorption peak 
appears in the frequency region between 350 cm~^ and 
550 cm~^. The frequency of its maximum increases with 
increasing doping; for optimally doped samples this fea- 
ture disappears. Secondly, at the same time as the peak 
forms, the infrared-active phonons in the frequency re- 
gion between 300 cm~^ and 700 cm~^ (in particular their 
strength and frequency) are strongly renormalized. This 
effect is most spectacular for the oxygen bond-bending 
mode at 320 cm~^, which involves the in- phase vibration 
of the plane oxygens against the Y-ion and the chain ions. 
For strongly underdoped YBa2Cu306.5 with Tc ^ 50 K, 
this mode loses most of its spectral weight and softens by 
almost 20 cm^^. Although the additional peak, and the 
related changes of the phonon peaks (phonon anomalies), 
start to develop above there is always a sharp increase 
of the peak magnitude below Tc |la] . Similar effects have 



also been reported for several other underdoped bilayer- 
compounds (see, e.g., Refs. |l^-p^ ]) and for hole-doped 
ladders in Sri4_3;CaaCu2404i |19||. 

Van der Marel et al. have suggested j2^ that the addi- 
tional peak around 450 cm~^ could be explained using a 
phenomenological model |^ of the dielectric response of 
superlattices with two superconducting layers (a bilayer) 
per unit cell. The model involves two kinds of Josephson 
junctions: inter-bilayer and intra-bilayer. As a conse- 
quence, the model dielectric function exhibits two zero 
crossings corresponding to two longitudinal plasmons: 
the inter-bilayer and the intra-bilayer one. In addition, 
it exhibits also a pole corresponding to a transverse op- 
tical plasmon. Van der Marel et al. pointed out that 
the additional peak in the spectra of underdoped YBCO 
may just correspond to the latter plasmon. Very recently, 
they have confirmed their suggestion by more quantita- 
tive considerations regarding the doping dependence of 
the peak position |Q. The details of the spectacular 
anomaly of the 320 cm~^ phonon mode, however, cannot 
be explained within the original form of their model. 

In the following we report a theoretical analysis of the 
additional peak and the phonon anomalies. We have ex- 
tended the model of van der Marel et al. by including 
the four phonons at 280, 320, 560, and 630 cm"! [||j23| 
in such a way that the extended model can account not 
only for the peak but also for the most striking phonon 
anomalies. The important new feature is that we take 
into account local electric fields acting on the ions par- 
ticipating in the above mentioned phonon modes. As we 
show below, the phonon anomalies are then simply due 
to dramatic changes of these local fields as the system 
becomes superconducting. 

Let us briefly introduce the model. The dielectric func- 
tion is written as 



e(w) = ei{uj) + ie2{uj) 
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where Soo is the interband dielectric function at fre- 
quencies somewhat above the phonon range, j„ are 
the induced currents, () means the volume average, 
and E is the average electric field along the c-axis. 
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The following currents have to be taken into account: 
the Josephson current between the planes of a bi- 
layer, ju — —iuJ£oXbi E^i, the Josephson current be- 
tween the bilayers, jint = -i(^£oXint Eint, the cur- 
rent due to the oxygen bending mode at 320 cm" ^, 



jp = —iujEoXpEiocP, and the current due to the other 
three infrared-active modes involving vibrations of ions 
located between the bilayers (apical oxygens and chain 

atoms), jA = -iujeoXaEiocA- Here 
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(3) 



are the susceptibilities that enter the model. The plasma 
frequencies of the intra-bilayer and the inter-bilayer 
Josephson plasmons are denoted as lou and uJint^ respec- 
tively. We do not attribute any physical interpretation to 
the Lorentzian terms in Eq. (2) that are designed solely to 
represent the featureless residual electronic background 
in the frequency range of interest (i.e., from 200 cm^^ 
to 700 cm^^) in a Kramers-Kronig consistent way. The 
response of the phonons is described by Lorentzian oscil- 
lators as usual. Further, Eu is the average electric field 
inside a bilayer, Emt is the average electric field between 
neighbouring bilayers, EiocP is the local field acting on 
the plane oxygens, and EiocA is the local field acting on 
the ions located between the bilayers. Note that by iden- 
tifying Eiocp with the field acting on the plane oxygens 
we have neglected the contributions of the other ions in- 
volved in the phonon (the Y-ion and the chain ions) . This 
seems to be a reasonable approximation, since the con- 
tribution of the plane o xygen s to the 320 cm~^ mode is 
known to be dominant |E2 23 1. 

The electric fields Ebi, Eint, EiocP, and EiocA can be 
obtained using the following set of equations: 



Ebi — E' + 



aXpEiocP 



£q£o 



Eint = E 



I PxpEiocP +1XaEi 



locA 



Eloc P — E' 



Elnr A = E' . 



-ILUK = Jint - ]bl 



E{dbi + dint) — Ebidbi + Eintdi 



(4) 

(5) 

(6) 

(7) 
(8) 
(9) 



containing two additional variables, k and E' . The for- 
mer represents the surface charge density of the copper- 
oxygen planes which alternates from one plane to the 



other whereas E' is the part of the average internal field 
E that is not due to the effects of k, xp and XA- The 
terms in Eqs. (4) and (6) containing k represent the fields 
generated by charge fluctuations between the planes. The 
terms in Eqs. (4) and (5) containing the phonon suscep- 
tibilities represent the fields generated by the displace- 
ments of the ions. The values of the numerical factors a, 
/3, and 7 (1.8, 0.8, 1.4) have been obtained using an elec- 
trostatical model While the feedback effects of the 
phonons on the electric fields have to be included in order 
to obtain the observed softening of the oxygen bending 
mode, they are not essential for explaining the spectral- 
weight anomalies. Equation (8) guarantees charge con- 
servation. The distances between the planes of a bilayer 
and between the neighbouring bilayers are denoted by du 
{dbi = 3.3 A) and dint (dint = 8.4 A), respectively. 

Figure 1 (a) shows the experimental spectra of the c- 
axis conductivity of YBa2Cu306.5 with Tc = 53K from 
Ref. p^. Figures 1(b), 1(c), and 1(d) show the data 
for (b) r = 300K, (c) T = 75K and (d) T = 4K 
together with the fits obtained by using the model ex- 
plained above. 

The values of the parameters used are summarized in 
Table 1. Those used in computing the room-temperature 
spectrum have been obtained by fitting the measured 
complex dielectric function from Ref. |jl^ (with ujbi — 0.0 
and ojint = 0.0). Those used in calculating the 4K spec- 
trum have also been obtained by fitting the data, except 
for £00, wp, and the oscillator strengths of the phonons 
(S'p, Si, 5*2, 5*3) which have been fixed at the room- 
temperature values. The appearance of the additional 
peak and the anomalies already at temperatures higher 
than Tc may be caused by pairing fluctuations within the 
bilayers. Motivated by this idea, we have fitted the 75K 
spectra in the same way as the 4K ones allowing only 
the upper plasma frequency (ujbi) to acquire a nonzero 
value. We shall comment on this point below. Note that 
the low-temperature value of ujint (220 cni^^) is rather 
close to the one obtained from reflectance measurements 
(204 cm~^ in Ref. ||l^) and that the screened value of luu 
(~ 500 cm^^) falls into the frequency region of a broad 
peak in the loss function (Im(l/e)) ||l^, which is a sig- 
nature of a longitudinal excitation. 
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FIG. 1. (a) Experimental spectra ol the c-axis conductivity, a = ai + i(J2, of YBa2Cu306.5 with Tc = 53K from Ref. [15]. 
Experimental data together with the fits obtained by using the present model for (b) T — 300 K, (c) T = 75 K, and (d) T — AK. 
The dotted lines represent the electronic contributions. 



TABLE I. Values of the parameters used in the present computation. The temperatures are given in K, the frequencies and 
the broadening parameters in cm^'^. 
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Si 
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72 
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6.5 


300 


5.15 








1280 


1000 


680 


250000 


1.3 


0.056 


0.087 


0.32 


390 


290 


553 


647 


21 


22 


29 


16 
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75 


5.15 
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450 


170 


810 


140000 


1.3 


0.056 


0.087 


0.32 


390 


285 


563 


655 
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14 


20 
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6.5 
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5.15 


1205 


220 


8.2 


2.7 


6000 


180000 


1.3 


0.056 


0.087 


0.32 


390 


284 


564 


656 


9 


12 


15 
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260 


5.15 
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220000 


1.3 


0.056 


0.200 


0.30 


390 


296 
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631 


13 


20 


15 
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6.8 
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5.15 


1779 


450 


13.8 


4.8 


2400 
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1.3 


0.056 


0.200 


0.30 


390 


288 


589 


646 


2.5 


16 


25 


20 


6.4 


250 


5.15 








310 


240 


860 


120000 


1.3 


0.060 


0.100 


0.42 


400 


286 


553 


650 


34 


24 


30 


15 


6.4 
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5.15 


949 


150 


560 


480 


520 


120000 


1.3 


0.060 


0.100 


0.42 


400 


280 


556 


654 


42 


24 


20 
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The values of the phonon frequencies are somewhat 
different from those which would result from a usual fit 
of the data (such as in Refs. §-|l|,||). This is because 
the susceptibilities of Eq. (3) represent response functions 
with respect to the local fields instead of the average field. 
In the absence of interlayer currents, the input frequen- 
cies would correspond to the LO-frequencies while the 
frequencies renormalized according to Eqs. (4)-(9) would 
correspond to the TO-ones. Our input frequency of the 
oxygen bending mode (390 cm~^) is close to the measured 
LO frequency |2|. 

It can be seen in Fig. 1 that the model is capable 
of providing a good fit of both the normal- and the 
superconducting-state data without changing the oscil- 
lator strengths of the phonons and without any change 
of the input frequency of the oxygen bending mode. It 
reproduces succesfully: (i) the appearance of the addi- 
tional peak, its position, broadening and magnitude; (ii) 
the loss of the spectral weight of the peak corresponding 



to the oxygen-bending mode and the pronounced soft- 
ening of this mode; (iii) the loss of the spectral weight 
of the peaks corresponding to the apical oxygen modes 
at 550 cm~^ and 630 cm^^ and the increase of their as- 
symetry. The intrinsic frequencies of the latter modes 
have to be slightly increased in order to reproduce the 
noticeable hardening of these modes. The dotted lines 
in Figures 1(b), 1(c) and 1(d) represent the results ob- 
tained after omitting the phonons in the fitted expre- 
sions {Sp = Si = S2 = S3 = 0.0). It appears that the 
plasmon peak collects the lost part of the normal-state 
spectral weight of the phonons. This, however, only ac- 
counts for a part of its spectral weight. In the absence 
of the phonons and for small values of the residual con- 
ductivities the spectral weight of the S — peak at a; = 
0.0 is Ss = {Tr/2)eo{dbi+ d^nt)(^uUjfnt / (dbii^Lt + d^ntUJu) 
and the spectral weight of the additional peak is 

Sp = (7r/2)£o(4irf»nt/(4i +d*»t))Ki -'^mt)V(4iwL + 

dint^bi)- For the values of the two plasma frequen- 
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cies given in Table 1 we obtain Ss = 1700 ^cm ^ and 
Sp = 10000 f]-icm-^. Note that both Ss and Sp belong 
to the spectral weight of the superconducting condensate. 
This should be taken into account in discussing the sum 
rules as, e.g., in Rcf. psj . 

Our model allows us to explain the decrease of the 



spectral weight of the phonons when going from the nor- 
mal to the superconducting state by using rather simple 
qualitative arguments. Neglecting the feedback effects of 
the phonons and the residual electronic background (rep- 
resented by the Lorentzians in Eq. (2)), the electric fields 
Ebi and are given by: 



Em = 



{dbl + dint)£int 17 77 _ {dbl + dint)ebl 



dbl^int + dint£bl 



dbl^int + dint^bl 



E, 



(10) 



where Ebi — £oo — ^u/^'^ ^^'^ ^« 



The low-temperature spectra of Ebi and Emt are shown 
as the solid lines in Fig. 2. In the frequency range of the 
oxygen bending mode, Ebi and Eint have opposite signs 
and similar magnitudes and the same holds for E^t and 
Ebi ■ As a consequence, the local field acting on the plane 
oxygens, which equals the average of the two fields Eint 
and Ebi (cf. Eqs. (4), (5), and (6)), can become rather 
small. The frequency range of the modes of the apical 
oxygen is close to the zero crossing of Ebi ■ Consequently, 
the local field acting on the apical oxygens, Eint, is rather 
small in this frequency region. It is the decrease of the 
local fields when going from the normal to the supercon- 
ducting state which is responsible for the spectral weight 
anomalies. The room- and low-temperature spectra of 
EiocP shown in Fig. 2 illustrate the above considerations. 




frequency (cm" ) 

FIG. 2. The approximate dielectric functions of the in- 
tra-bilayer and the inter-bilayer regions, eti and ei„t defined in 
the text (solid lines) . The room- and low-temperature spectra 
of the local field EiocP acting on the plane oxygens (dashed 
and dotted lines). 



Our model is also capable of explaining the experi- 
mentally observed doping dependence of the additional 
peak and the anomaly of the oxygen bending mode. As 
the doping increases the peak shifts towards higher fre- 
quencies and it becomes broader and less pronounced 
(see Fig. 10 of Ref. and Fig. 3 of Ref. |l|l). Both 



these trends can be easily understood. The first one is 
due to the progressive increase of the plasma frequencies 
with hole doping, which is related both to the increase 
of the condensate density and to the reduction of the 
charge confinement |14|. The second one is due to the 
fact that the broadening is proportional to the residual 
backround conductivities which increase with increasing 
doping. In addition, the size of the spectral gap {2l^max) 
decreases with hole doping jl4| and eventually falls below 
the energy of the transverse optical plasmon around opti- 
mum doping. The phonon anomaly appears in the same 
range of doping as the additional peak. For y around 6.8, 
the spectral weight from the high-frequency side of the 
phonon peak moves into the additional peak (see Figs. 10 

(a) and (b) of Ref. |T^) as the temperature is lowered. 
For y around 6.6 we find the most pronounced anomaly 
(see the experimental data of Fig. 1). For even lower 
doping levels the additional peak and the phonon merge 
together forming a single highly-assymetric structure (see 
Fig. 3 (d) of Ref. [|llj). These trends can be understood 
using arguments similar to those presented above (see the 
discussion related to Fig. 2) and can be well reproduced 
using the model. This is demonstrated in Fig. 3 which 
displays the experimental spectra of the c-axis conduc- 
tivity of YBaaCugOy with (a) y = 6.4 (T^ = 25 K) and 

(b) y = 6.8 (Tc = 80 K) from Ref. ||l|l together with the 
fits (Fig. 3(c) and 3(d), respectively). 

We discuss next the peculiar temperature dependence 
of the additional peak and the phonon anomalies. The 
proximity of the onset temperature of the anomalies and 
the onset temperature of the spin gap (T*) observed in 
nuclear magnetic resonance experiments has provoked 
several speculations [|l^,0 that the anomalies are due 
to the coupling of the phonons to spin excitations. From 
the fact that we are able to fit the data for temperatures 
between Tc and T* (see Fig. 1(c)) we infer that in this 
temperature-range the intra-bilayer plasmon is already 
developed. This suggests that many of the electronic and 
possibly also structural anomalies starting below T* (see, 
e.g., Ref. j2^) are caused by pairing fluctuations within 
the bilayers. 

In summary, we have extended the phenomenological 
model of Van der Marel et al. involving inter-bilayer and 
intra-bilayer Josephson junctions by including phonons 
and local field effects. 
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FIG. 3. Experimental spectra of the c-axis conductivity 
of YBaaCusOj; with (a) y = 6.4 (T^ = 25 K) and (b) y = 6.8 
[Tc = 80 K) from Ref. [15]. (c) and (d) Fits of these spectra 
obtained by using the present model. 



The model allows us to explain not only the additional 
broad peak around 450 cm~^ but also the spectacular 
anomaly of the oxygen bending mode at 320 cm~^ and 
the spectral weight anomalies of the apical oxygen modes 
at 550 and 630 cm^^. Our results indicate that also the 
closely spaced copper-oxygen planes of underdoped bi- 
layer cuprates are weakly (Josephson) coupled. These 
findings provide support for the conjecture Q that the 
c-axis dynamics of the cuprates (at least the underdoped 
ones) is dictated by the unconventional properties of the 
ground state of the electronic system of the planes. We 
suggest that the onset of the anomalies at temperatures 
significantly higher than Tc may be caused by pairing 
fluctuations within the bilayers. 
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